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The intrinsic tau neutrino flux from cosmological and astrophysical sources has usually been considered negligible in com- 
parison to the electron and muon neutrino fluxes. However, the inclusion of the tau neutrino component coming from hadronic 
decay at the source can significantly modify the tau neutrino spectrum expected at Earth. We report our results on the high 
energy tau neutrino production and its implications for the observation of high energy neutrino events. 



1. INTRODUCTION 

Presently much effort is being expended to resolve 
the flavor of the incoming neutrino flux (for a review 
see [|ljj). One motivation for doing this is that the 
absolute sensitivity, or relative sensitivity compared 
with background, may be significantly greater for one 
neutrino species than for other species Another 
primary motivation for identifying the neutrino fla- 
vor is to test the hypothesis that neutrinos undergo 
flavor oscillation as they propagate from source to 
Earth. Because the non-oscillation production of v T 
is assumed to be negligible relative to v e and v „ pro- 
duction, and the cascade v T produced by collisions 
of primary v e and flux with the relic cosmic neu- 
trino background is orders of magnitude less than the 
primary v e and flux, any detected v T component 
from a source outside the Solar System is currently 
expected to be indicative of oscillation. In this paper, 
we briefly report some results of our investigation into 
high energy v T production, presented elsewhere [^J. 
We illustrate that the v T component at the high energy 
end in the hadronic decays is significantly higher than 
previously assumed and that this can have observa- 
tional consequences. 

2. NEUTRINO PRODUCTION 

High energy neutrinos are expected to be produced 
by several cosmological and astrophysical sources, 
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among others: active galactic nuclei (AGN) |Q|; topo- 
logical defects (TD) such as superconducting, or- 
dinary or VHS cosmic strings [§-0]; supermassive 
gauge and scalar particle (X-particle) decay or anni- 
hilation and Hawking evaporation of primordial 
black holes (PBH) [pLpHj. Also, neutrinos can origi- 
nate from the decay of photoproduced hadrons on the 
cosmic background radiation (CMB). In all these sce- 
narios the production of tau neutrinos has been as- 
sumed to be negligible compared to the production of 
electron and muon neutrinos. 

We focus our attention on those cases where, irre- 
spective of the scenario, there are sufficiently ener- 
getic interactions that quarks and gluons which frag- 
ment into jets of hadrons are expected to be produced. 
The spectra of these hadrons should be very similar 
to those measured in e + e~ — > qq — > hadrons events 
in colliders. As a consequence, the final state par- 
ticle distributions in these scenarios are expected to 
be dominated by hadronic decays at the source. In 
these decays over 90% of the final state products will 
be 7T°, 7T + , and ir~, with the remainder mainly nu- 
cleons which decay into protons and antiprotons. On 
astrophysical timescales, the 7r ± decay into v e and 
and e ± . The final cluster states (pions and nucleons) 
in QCD jets at accelerator energies can be loosely ap- 
proximated by the Hill, Schramm, and Walker (HS W) 
fragmentation function 

f ->-. 

where x = E/mj and mj is the total energy of the 
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decaying jet [jg]. This distribution continues down to 
E ~ 1 GeV. When convolved with the 7r ± decay, Eq. 
([[]) leads to a similarly dominant _E~ 3 / 2 term in the 
j/ e and spectra at x < 0.1 and dN Ueill /dE — > as 
a; — » 1. (For a full derivation of the v e and spectra 
using fragmentation function ([[]) see [Q].) 

Because zy T production is suppressed compared 
with other species, it has been assumed in previous 
astrophysical and cosmological flux calculations that 
the v T spectrum from hadronic jets is orders of magni- 
tude less than the v e and Vu spectra at all x. This is not 
so, once the energy of the jets surpasses the tau lep- 
ton and heavy quark masses. While indeed the total 
number of v T produced per jet is less than of the 
total number of v e and v^, the high x tau neutrinos 
are predominantly produced by the initial decays of 
the heavier quarks with shorter lifetimes (the greatest 
contribution comes from the t quark decay) and the v e 
and are produced by the final state cluster decays 
of the much lighter pions. This leads to significantly 
greater relative contribution from v T at high x than 
previously assumed. The fragmentation distribution 
(0) is no longer relevant for the tau neutrino. 

In Fig. 1 we show the v T spectrum, together with 
the v e and spectra, generated by the decay of 10 
TeV qq jets. To simulate these spectra we used the 
QCD event generator HERWIG [|l2|| and the process 
e + e~ — > qiqi (g), i = all q flavors. Consistent spectra 
are obtained with PYTHIA/JETSET (see [|]). Note 
that in the region 0.1 < x < 1, the tau neutrinos make 
up more than one tenth of the total neutrino contribu- 
tion. Note also that below x < 0.1, dN Ur /dE falls 
off with roughly an E^ 1 ! 2 slope, and not the i?~ 3 / 2 
slope of the v e and spectra. From our 300 GeV 
- 75 TeV simulations, we find that the v T spectrum 
generated by qq jet decay can be parametrized as 
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per jet. 

Here we apply our results to two scenarios (for 
an expanded treatment see Q]). The first one is the 
VHS cosmic string scenario. In the VHS scenario, 1 
GeV - UHE particle fluxes are generated by the de- 
cay of supermassive scalar and gauge particles emit- 
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Figure 1. The neutrino spectra, dN/dE vs x = 
E/mj, generated by the decay of raj = 10 TeV qq 
jets. The dashed line represents v^, the dotted line 
represents v e , and the solid line represents v r . 



ted by the long strings over the age of the Universe. 
In Fig. 2 we show the neutrino spectra expected at 
Earth from VHS strings with a mass per unit length 
of G/i = 10~ 8 [0]. The partial v T component which 
is solely produced in the collisions of the primary v e 
and with the relic cosmic 1.9°X neutrino back- 
ground [ jl3t l is also shown in Fig. 2. This collision- 
produced v T component is the only one presented in 
previous cosmic string papers. As can be seen, our 
results give a significant increase in the expected v T 
signal at the highest energies. In the second scenario, 
the evaporation of PBHs, we calculate the neutrino 
flux from the Hawking evaporation of a Tbh ~ 100 
GeV PBH using HERWIG and following the method 
of Ref. [[Top. The instantaneous emission per degree 
of freedom prior to decay is given by 

-i-lf 
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dQdt 

where T s is the absorption probability for species, s 
is the spin, and Tbh is the black hole temperature. 
The neutrino flux from the evaporating PBH including 
hadronic decays is shown in Fig. 3. 

3. CONCLUSIONS 

Previous work has assumed the v T component in 
TeV - UHE hadronic decays to be negligible com- 
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Figure 2. The neutrino spectra from VHS cosmic 
strings with a mass per unit length of G[i — 10 -8 . 
The dash-dotted curve is the v T component produced 
by cascades off the relic neutrino background [|l3[] and 
the thick solid line is the the v T component produced 
by hadronic decays. 



pared with the v e and Vn components, at all neu- 
trino energies. We find that this is not so once the 
decays are sufficiently energetic to include the heav- 
ier quarks. In particular, for neutrino energies in the 
decade below the energy of the initial decaying par- 
ticle, the tau neutrino component is of similar mag- 
nitude to the v e and components. Below these 
energies the v T spectrum exhibits a slope of slightly 
less than E~ x l 2 compared with the E~ 3 / 2 slope of 
the v e and spectra. This analysis modifies the ex- 
pected spectra in many astrophysical and cosmologi- 
cal high energy neutrino production scenarios. As a 
consequence, the observation of a significant v T to 
ratio at a given energy in high energy neutrino tele- 
scopes and detectors may be due to hadronic decay at 
the source and not — > v T oscillation in transit. 
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